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SYMMETRIES IN LAMINATED 
COMPOSITE PLATES 

AHMED K. NOOR 

George Wash4ftgton University Center 
at 

Langley fteseerch Center 


Abstract 


Ire ail’ferer.t types of sytanetry exhibited by laminates anisctrttic 
fibrous composite plates are identified and contrasted vith the symmetries 
of isotropic and homogeneous crthotrcpic plates. The effects cf varia- 
'.ions in the fiber orientation and the stacking sequence cf the layers on 
the symmetries exhibited by composite plates are discussed. Both the 
linear and geometrically nonlinear responses of the plates are considered. 

A simple procedure is presented for exploiting the symmetries in the finite 
element analysis. Examples are given of square, skew and polygonal plates 
where ise of symmetry concepts can significantly reduce the scope anu cost 
cf tr.alysis. 


Introduction 


tn the analysis of isotropic and orthotropic plates, ccsputaticns can 
be substantially reduced in scope and cost if certain symmetries exist in 
geometry, material properties, loading and boundary conditions. For 
-im inated composite plates cade of anisotropic materials s;cocetry proper- 
ties of a different nature exist and except for a recent study by the 
author [Ref. l], these symmetries have not been utilized in finite element 
ana.yses. The present study focuses on this problem and is an extension 
cf t.'.e work reported in Ref. 1. More specifically, the cb.;ectives of this 
paper are: a) to identify the different types of symmetry exhibited by 

the cocEcnly-used fibrous composite plates in both the linear and non- 
linear regimes and b) to present a simple procedure for exploiting these 
symmetries in the finite element analysis. 

Tne plate is assumed to consist of a numoer cf layers bonded together, 
lacn luye" has arbitrary thickness, elastic properties and fiber orlenta- 
•.I’.n vith respect to the plate axes. The analytical formuiaticr. is based 
n a ge-.metricaliy nonlinear von-Karman type plate theory with the effects 
:f traT.s/erse shear doformation, .misotropie material behavior anu 
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b«ndin?-exter.s ior.il coupling included. A displacener.t finite elenent 
fomulation is assuced with the findanental unknowns consisting of the 
displaceaent and rotation components of the middxe p.ane of the plate 
(ur.. w and p.,). For convenience, the nodal parameters are taken to be 
the five generalized displacements at each node. The sign convention for 
the different plate quantities is shown in Fig. 1. It shoolc be enp.ha- 
sized, however, that the proced-ure outlined herein for exploiting the 
syiaetries in the finite elenent analysis can be adapted for use wiU. ^ 
other plate theories and other finite element models (e.g. mixes noueis.. 

Examples are given of sq'uare, skew and polygonal cosposite plates 
where use of sj-=etr>- concepts can significantly reduce tae conputational 
cost of the analysis. 


■■unianental definitions of symmetry, symmetr;.- e.eaen.s 
atiins are reviewed in Ref. 2. The "axiom of sysattry 
pef, 2 when applied to composite plates reads; .ji 
tlate exhibiting certain types of symmetr;.- and a syste 
hits the same types of symmetry as those of the p_ate, 
htained will exhibit the sane types of sytmoetry as the 
.n the plate and the loading system. Here the symmetry 
srs to the sjT=etry cf a) plate geometry, h' .^.inatic 
:. fiber orientation and stacking of layers) ano c' b: 


of the plate is 
s the distance 


a ii/mmctti/ t'uxmioWdtict: 
:s the iynrctly sf t 


sjTsmetry frequently ( 
bilateral) symmetry; 
,ry. The characteris' 


e plate exhibits reflection symmetry \ 
try) with respect to a plane, a x^ ^ 
ougr.t into an equivalent configuratit: 
Obviously, the loading on the plate i 
possess mirror symmetry with respect 
E which are skev-symnetrically lamina" 
a.ne [.=.efs. 1 a.-.a 5) often exl.itiu a d 


the boundary 
the same plane. 

. with respect to 
'erent type cf 


in that plane 
conditions mu 
Composite p_a 
their miiile 
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- s -.arelv. reflection sj-cnetry with respect to a 

• s,-*^ — V ^loTA rhus ze tr^kVr 


cilater^ s:-=e ry be into 

a X'_ + 0 X2 , If. ...e*. _ ,ar>c ,„t„tier. abo-t tne line. 


« i,^iva^?nt ccnfS^^a^I- th^:;u^h'a Iso^ rotation abo-t tbe line, 
a X- ^ X''* 


* ~ve tvc tvnes of sy^etry, reflection in a plane ^i^reflecticn^in a 
— b t» called reflection sytsaetries c- .-«, = 
r:i:/:r,e'' "rs^ectively. The s:/=etry relations for tne ;*•«,; -=P°- 

;:i!^:., reflection t:.Te 1, is the better "-^^t-e-ric 

■■,iv 4 ."h V... ^so^rocic and honcgenecus orthotrcpic pxa..o *“* "t 
;«:e:;Ssr ;i";%hat the sy»etry relations^ " uT -d 

^.e.. in-plane loading exponents typ^s I 

:hr”h;r^d:“th;‘^fcx.etry relations for the benoin._^oup 

.‘-.e., transverse *cal*..(^ p , _ ♦••ar.sverse s’r.ear 

conpenents ^,r,er In sign fren the 

stress resu-.ar..s -a- I*; s^^^ed differently, the 

for’ the bending group in reflection type II are 
identical vith the antisytaetry relations in reflection ./?e 


Composite plates can possess one r'^ a* e^ with'^ infinite 

tion sictmetry. Circular plates plates 

number of planes of reflection si-mmetry. '1 ir. 

exhibiting reflection syraitetries cf tne -ypes 
figures 3a and 3b. ' 


■'gures sa and 3b show conto-ur plots for the sut^c^er 

Ii.pi;:!ment components -“^r/rig^^C’a^ 

to -uniform transverse respe^velv . :iote that the 


= Xj). Cr. the other nano. ^ ,_5-5 rctaticn about 

can be bro-ugnt into co*nci-er.ce ^ g ^ equivalent to a 

the line cf reflection symmetry (line *1 ‘ *2J- ^=“ ^ 

-rror reflection in the plane xj - ^ through l8CO about 


mirror reflection in tne p-ane xi - xg* 'tv^c'ugh iBC® about 

mirror reflection in plar.e « x or a ro-aticn t xc-gn x 

line X-, . x„ leave: the contour plots for w. in fig-es sa ar - . 


unchanged. 


rot.atlcr.al sv-ixietry 


A conpo.lt,, pl.t, 1 . ••IJ t- «l.iwt "t.tlon^ 
to „ „i."»rn.i to It, n»' ‘"“'“i; =; St:;i“ 

:riL^j;n!"rori s:ri™ 

Z‘i'iS,‘t“'wr‘tl». n..o,l,t,i .lU, It, 
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aanely rota-sisns of 2*:, it, 2^/3, ... 2it/n (see Ref. 2/. The s;.'=r,etr>' 
relations for the loading cocponents, generalised dispiaceisentc. ar.d stret.. 
resultants for rotational sj’tsaetry vith n»2 and n«» are giver, in tables 1 
and II. The axis of rotation is assumed te coincide vith the Xj-axis. 

Note that for n»l the syrscetri' relations for the tendlrj* ctreo. resul- 
tants in ouaoi-isotrsf ic cor.posite plates f see Ref. - differ -r. ai.'r. 
from the corresponding ones fer isotropic plates. 


An important special ease of rotational symmetry is *he cere :f 
complete or axial s;.-smetry (r.=*); which is exhibiteo t;.- Cw^mpoeite circu- 
lar plates whose elastic characteristics and boundar;.- conditions are 
independent of the circumferential coordinate (i.e. axia;nnr.etric . The 
response of the plate will also be axisjnsmetric . Kewever, in contra- 
distinction to isotropic and orthotropic plates, the two sets of dis- 
placements and stress resultants (uj, dij, Tr,„, Mu, ^,)*and (u. , w, C- , 

, , S_., M,, M-, 1- ), though axisymmetriSt ar^^coupleu, i.el tney 
d6‘‘’nct^all ^Sniss^under an;,’ cf the axisymnetric loadings shown in Fig. 2 
(see Ref. 6). 


It should be mentionec that reflection s;.nnmetr;.’ type II car. be con- 
sidered as a rotational symmetry, n*2, with respect to a line in the 
middle plane of the plate. 


Inversicr. s~.~uov*-r:.’ 


- A composite plate is said to exhibit inversion s^.’metr;.’ ’.’ith respec*. 
to an axis x- normal to its plane if It car. be Drought into a.', equ. valent 
coniiguratlon through 18C° rotation about the axis. This amounts tc 
changing the coordinates x^ of each material point cf the p_ate into 
-X . Tr.c x--axis is called the axis cf symmetr;.- and its intersection 
wi^. the miille plane of the plate is caliea the center of st’mmetr.v . In- 
version symmetry ear. be thought of as a reflection through a point, vio. 
the center of symmetry; or a rotational symmetry with r.=T (i.t. rc-.ation 
angle*") . 


Tr.e center of st’mmetry is located at the intersection cf — nes of 
geometric s;.c:cre-r;’ or. the middle plane of t.ne plate. For p-ateo wotn sx* 
plar.fcrr. the center -of symmetry is located at the intersection of tae f.’: 
diagonals. An example of a cccpxjsite skew plate sublectcd to -r.iform 
loading and exhibdtir.g inversion type symmetry is given in Fig. -. Tne 
inversion symmetry relations listed in tab^e II are clearly oevious in 
the contour plots of Fig. 1. 

It may be mer.’ijr.ed that for certain geometries arid- or loai-ngt t;.; 
response of isotropic (and orthotropic) plates exnioits Ir.versi.:. syt.- 
netry simi-ar to tr.at exnibitea by anisotropic vcoeposite, p.it-u. 
Examples of tnese si'uatior.s include plates witr. skew rianJ:."m and plate: 
subjected to edge shear. 
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ionllr.ear 


Th« ar^inento of reflection s>Tn:etry in plane (re: 
ersior. ana axial s^Tnetry apply tc both the linear i 
linear reBponses of the plate. Cr. the ether hat., 
h r.=- for quasi-lsotropic plates and reflection s;.r; 
flection type II ' holi only for the linear respcr.se 
lalnec. for th» latter sjTasetry by the fact that in ■ 
linear case, a conpcslte plate behaves like a s.nall 
at ion of 180° about the axis of reflection sytaretr;' 
face) does not bring the shell into an equivalent c 


For anti jjTjnetrlc (or skev-syemetrie ) .oadings 
ate will be antlsynmetric (or skew-syiuaetr Ic ) . r 
on saev-syncetry the transverse displacement w, 
sultar.ts and the Lendlr.g stress resultants 

r.ter of symSitry. 


of the 
inver- 

• 638 “ 

. the 


Any loading system can be decomposed into synccetric and antisyTmnetrlc 
(or skew-symmetric) components. For linear problens, significant computa- 
tional advantages result from such a decomposition. 


In free vibration problems of composite plates, the mode shapes car. 
be obtained by using a portion of the plate and applying the frir dif- 
ferent combinations of symmetry and antioynr.etry (cr 3 /.eu-symnetry . cor.di 
cions along the internal boundariee. Such an apprcach was found to 
significantly reduce the computational effort in the case of isotropic 
and orthocropio plates. 


Composite plates having the sane symmetric geometry, loading and 
tcu.ndary conditions, but llfferert laminations car. exhibit different types 
of symmetry. The lamination para, cters which have t.ae strongest effect on 
symmetries are fiber orientation of Individ'ial layers a.nd stackir.g sequence 
of different layers. To llluatrate the effect of tr.ese two parameters on 
the sjTsmetries, table III gives Uie oyTmetry trar.sfcrmatlonr for conposite 
square plates having t.he following laminations (fc~ a deacription of the 
characteristics of these laminates sec Pef. •*). 
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ii) Midplane antlsj'sr.etrlc lar.lnatea (cross-ply, or.<;le-pl: 
and angle-ply with 


ill) Quas 
9»+1*5/0/90/-1*5 


sotropic, four-layered lar.inate with 


The loading on the plate and the boundary conditions are ass-jned to be 
sjxasetrle. For the purpose of ocnparison, the cioe cf ‘he roithen-iticu- 
model required for the ana-ysis as well as the siT.-r.etr;- trtr.tftrruitior.o 
for Isotropic and homogeneous orthotroplc square plates are also gi-/en in 
table III. The SiTsetr;.' relaticns for the generalised displacements and 
stress resultants, implied by t.he Si-mmetry transformations of table III, 
are given In table II. 


Figures 5 through ' show contour plots for stress resultar.ts and 
generalized displacements in composite square plates with symmetric, 
antisymmetric and quasi-isotrcpic laminations. The plates have clamped 
edges and are subjected to uniform transverse loading. Tne symmetry 
relations listed in table II are clearly obvious in these flrures. h'ct 
that for quasi-isotroplc plates (Fig. 7) there are no planes or iinec o 
reflection Sittmetry and the symmetry gro-up of t.he p.ate inc.uaes rctati* 
through 90° (and 1:0°) abc-..t the Xj-autis. If these s;-mmetries are 
utilized, only one quadrant of the plate needs to be ar.alyred. 


reflection in 
response of * 
response. Th 
rotational sy 


hosen wuj( exr.it: 
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s;"=»‘.ry as thcs* of the plate. This Dear.s that the ear. te brought 
Ir.tt self-''oinsIler.ce bj each of the symmetry trarisf:matIor.s .;sei for 
the plate (rotation about axis of revolution; nirror reflection in each of 
t.oe planes or lines of reflection syrsmetry ar.i inversion through each of 
*he centers of s;.Tme";ry). 


Ir.ietendent no-ies ar.i elenentc 

After the finite element grid is selected, the next step is to iden- 
*ify the nodes assccitted with the independent degre*.: -f freedom of the 
plate. Ihese noses will .henceforth be referred to as ■ijidepe.rJijit ncdzi • 
rne t .r nodes in the plate will be called ie,oetufc«* t’.odu • -he depend- 
ent r des will be designated by the symbols, S, I, P :r L according to 
v.neti.er they are obtained from inuependent nodes by rotation about axis 
;f revolution, inversion through center of symmetry, reflection in a 
plane or reflection in a line, respectively. 

The minimum number of finite elements which, by successive appllca- 
ions of the siTtmetry transformations can cover the v.-.cle plate will be 
rcCerrea to as the set a'-oidet-'endctit c£e/«eiUi. Other elements in the 
p.ate will be called dependent elements. According t: this definition an 
iniependent element caru-.c* be obtained from other independent t.emei.ts by 
symmetry transformations (e.g. rotation about axis of revolution, reflec- 
tion or inversion). 

The mwi! tiptccity of or. element is defined as the r.-omber of times this 
e.ement con appear in the plate by symmetry transformations. It is equal 
t; one plus the number of dependent elements that car. be obtained from 
that Independent element by sjTsmetry transformations. Tne size of the 
finite element model ani'the number of simultaneous algebraic equations 
required in the analysis are governed by the number of l.ndependent ele- 
ments end nodes, respectively. Figures 9 and 1C show tr.e independe.at 
nodes and elements for plates with reflection and inversion symmetries 
modelled by rectangular and parallelogram elements. :;ote that the 
Inoependent nodes and elements are not unique. y given and con be selected 
in many different ways (see Fig. 10). 


Procedure for exploiting the symmetry 


A simple procedure for exploiting the sytssetry in the finite element 
analysis of composite plates is outlined in Fig. 11. Tne key elements In 
this procedure are discussed subsequently. 


The following modifications are made to the irput data: 

1_. The total number of nodes and elements are set equal to the 
number of independent nodes and elements, resieotlvely. As 
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a constquance of this, tht coordinatat of the Independent ■ 
nodes only are input. Also, the connectivity lists of the g 
independent elenents only are prescribed. 

a 

£. The aultlplielty of each eleaent aust be prescribed. This 
can be input as an additional entry in the eleaent con- 
nectivity table. 

a 

- — If any of the independent elesents has nodes of the type_R, 

1 , P or L , they oust be identified, in the eleae nt eon-^ — 
nectlvity list (see Ref. 1). 


Eleaent r.tdu:.eB 


The eleaental matrices for the Independent elenents c.nly are gen- 
erated. For elenents with nodes lying or. an axis of reflection synaet/y 
a e b X 2 (a, b are nonzero constants) or vith dependent nodes of the 
t}-pe R, P or L , the stiffness, nasa anif load matrices are modified 
by simple transformations of the form: 



O (-) 


[K] - ir]*[K] tr] y 
[)J] • [rn«] m 

s^[P] - _ I 


where ’K], [ttl ard iP] are the element stiffness, mass ana load 
matrices, tT] is a transformation matrix and a car over a matrix deletes 
trar.sfcrmed (or modified) matrix. O c ■ 


If the nodal partr.eters are llstea node by aede, the matrix [F’ will 
have a block diagciial structure, i.e. 



where m is the r.vimber of nodes in the element. For shear-flexible 
ele-er.*s with ncsal -arameters consisting cf the olsrlacenent and rotaticr. 
cocj-.;.er.*i, the .7], are 5x5 subcatrlces . The itrmi ^ these 

subma'rlres at I , Tar.l L nodes are giver. In Atperxis A. 

Finally, the elemental matrices (stiffness, mass ar.. ..taa matrices, 
cf each if tnu independent elements are multiplied by the T..ltipllcity cf 
that element. This step can be bypassed if all '.te ihdepen.ent eleser.it 
have the sac- multiplicity. ~ 






IF*' 




I lij 
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The eleKen*.a- Mtrlees for the inaepender.t eleser.ts only »re 

Zr. the asjer.tly proceie no distinction Is rede between the 
ieper-ient nodes .e.fi., S , I , P end L ) within the independent elesents 
.01 thslr corresponding Independent nodes. This can te aecoaplished by 
seating the r.-.r.ters cf the dependent nodes equal to the nuskers ef the 
c:r re spending independent nodes prior to the assembly and then the sssem- 


the usua> r.anner. 


stion of 


The boundary scr.ditions along the edges cf independent elements are 
applied in the usual sar.ner. In addition, the sjTBLetry and sitcw-sj-mcetry 
conditions muil be applied at each of the centers of Inversion s;.-=cetry and 
alcng the axes and lines of reflection syssaetr'- Th.e s:.Tsnetry and skew- 
?;’msetry conditions for snear flexible elements are listed in table IV. 

The resulting reduced set of algebraic equations is solved for the inde- 
pendent degrees cf freeocn. •• 


The foregoing procedure for exploiting the symmetr;/ in the finite 
element analysis of cempesite plates has been implementei in a small 
research-oriented finite element system and was found to result in cen- 
lilersble savings in the compute time required for the analysis. The 
savings ore primarily due to: a) the generation tnd assembly of the 

elemental matrices of Independent elements only and b) the reduction in 
the sice of the resulting system of simultaneous algetrsic equations 
Ithelr number, and possibly, their bandwidth). The savings can be 
particularly significant for eigenvalue and nonlinear problems. 


Concluding Remarks 


The different types cf symmetry commonly exhibited by composite 
plates for various loadings and boundary conditions, are identified and 
contrasted with tne sitaetrles of isotropic and orthctrcpic plates. Both 
the linear and gecnetrlcaily nonlinear responses of the plate are ecn- 
sidered. A simple procedure is presented for exploiting the s>caietries 
.n finite element analysis. The analytical fcmulatlcr. is based tn a 
geometrically nonlinear plate theory with the effects ef material 
arlsctrcpy, ber.dlng-extenaional coupling and tranaverse shear deformation 
inel..;ed. A iisplaeer.er.t (stiffneaa) finite element form-ulatlon is 
ttaumed wutn tne nodal parameters eonaistlng cf the disp.acement and 
ritetl&n components of tl.e plate middle plane. However, the procedure 
.utllned for exploiting the eynmetry can be readily used with ether plate 
tr.egrlee and other finite element models. 


.y-rplei are given cf cespoaile square, skew and polygonal plates 
-iiere use •.; car. slgr.ificai.liy reduce the n.j:ler uf Independent 

-•/r-jis f frei-lcr., ur.d \ -,ce the computational time repulred for their 
fir it* fiiaet.t i.p : ; . 
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of lit. '.Ai) and iA2) are iiitlpllei by a m-nit t.RR. 












AppcndU B - :iot«Xicn 

fcllcvlnc tj-Bbcli »r« u*m in this paper: 
a. ,a^ ai<*e len^hs of the plate 

h thichr.asi of the plate 

~.r!«ir.al and codified stlffneii catrlcee of the plate elexen' 
c r.'Aber of r.odei in the elenent 

[uJiLm] original and codified case natrlcai of the plate elesent 

M^g ber.iir4! streie reeultonte 

5 j extetielonal (in-plone) etreee reiultonti 

P„tP external load inteneitiee in the x^ and x, direction 

c ad 

original and codified load catricea of the plate elecent 

transverse shear stress resultants 

a 

J^,w displaceser.t oceponents in the coordinate directlcns 

u^,a_ lisrlac«ent cocponents alonr and norr.al to the axis of 

reflection s;.T3r.etry 

X ,x, Cartesian coordinate systen (x, ncnsal to the siddle plane 

* of the date) 

rotation cocponer.ts 

K,K tracing constants which con have the valves *1 and -1 

6 fiber orientation ar.gle 

[r] transfcmation satrix 

[r], $x^ '.rant format ion subcatrices 

I ] rectangular or siiuore catrix 

H,J,P and i denote dependent nodee obtained froe toe independent nodee 
by rotat.on. Inversion reflection ir a plan* cr reflee- 
tiou in a line. 

The ror.ge uf the Creek indices a,B is 1,^. 
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original 
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Table IV 


Syainet.ry and antisyrae^ry I r skew-syianetry) conditions for shear-flexible 
displaceiLent finite eleoent codeia. 
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b) Fiber onenlafion 0/90 


Figure 3 Contour ploU for ditplocemenft. T«o-loyered compoeite plate 
with one plane (or line) of reflection eymmetry tubjected to 
uniform trontaerte loading 
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a) Cross ply 0/90 


Contour plots for displocsmtnts. Two* loytred ontisymmstrlcdlly 
lominotsd clamptd squort plots subiscisd to uniform 
tronsvsrst loading. 
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Fibor oriintotion 4S/0/90A4S 


Figurt 7 Contour plot* (or diiploctmtnt* and (trttd rttultonti Four 
laytrdd au»i - itotropic SQuort plo'« witn (ibtr onenlotion 
45/0/90/- 45 *ub|«cl*d to uniform trontvorst looding. 
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Plates with two planes lor lines) of reflection symmetry x, = 0. x, = a Symmetry types I and la. 


ty Plates with two planes tor lines) of reflection symmetry x = ♦ x_ Symmetry type 1 1 


c) Plates with four planes lor lines) of reflection symmetry x. = 0. x. 
Symmetry types IV and I Va. * 


Figure S Examples ol plates with two and four planes lor linesi of reflection symmetry ana 
their possible symmetry types isee table III). 
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Modify input data to account for 

•Dependent nodes and elements 
•Symmetries in the model 


Elemental i routines 


Generate elemental matrices for independent 
elements only 

Modify elemental matrices to account for 
symmetries 


Assembly! routines 


Assemble elemental matrices for independent 
elements only 

Dependent nodes in independent elements 
are treated as independent nodes 


Boundary and f co nstraint ^n^iti^n^ 


Apply symmetry conditions at 
•Centers ol symmetry 

•Nodes lying on axes of reflection symmetry 


Solution! m odule 


Solve ’’educed set o* ‘K:'j‘*’Ons 


Figure 11 : Flow chart of the prccedu'p for exp oiling ^.mi 
finite element analysis. 




